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Abstract 

A Coherent Potential Approximation is developed for s-wave and d-wave 
superconductivity in disordered systems. We show that the CPA formalism 
reproduces the standard pair-breaking formula, the self-consistent Born Ap- 
proximation and the self-consistent T-matrix approximation in the appropri- 
ate limits. We implement the theory and compute T c for s-wave and d-wave 
pairing using an attractive nearest neighbor Hubbard model featuring both 
binary alloy disorder and a uniform distribution of scattering site potentials. 
We determine the density of states and examine its consequences for low tem- 
perature heat capacity. We find that our results are in qualitative agreement 



with measurements on Zn doped YBCO superconductors. 
Pacs. 74.62.Dh, 74.20-z 
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I. INTRODUCTION 



A treatment of disorder is an essential part of the theory of superconductivity. After 
all, one must explain why impurity scattering does not cause resistance. Thus it is natural 
that as evidence for novel superconducting states multiplies the foundations of the subject, 
due mainly to Anderson^ and Abrikosov and Gorkovi'i, are being re-examined. The ex- 
periments which stimulate most strongly the current revival of interest in the problem are 
those on the high temperature superconductors!, which are now universally regarded as 
'c?-wave superconductors'!, and those involving some of the heavy fermion systems which 
display signs of 'p-wave' pairingi. In what follows we wish to contribute to the theoretical 
discussioniHli of the issues raised by these very interesting developments. 

The case of classic, 's-wave', superconductors is by now well understood. If the per- 
turbation does not break time reversal symmetry and the coherence length is sufficiently 
long, so that the pairing potential A does not fluctuate, the Anderson TheoremS guarantees 
that there is an absolute gap in the quasi-particle spectrum and the main effect of disor- 
der is that the density of normal states in the gap equation is replaced by its average over 
configurations^ On the other hand if the perturbation breaks time reversal invariance, as 
is the case with paramagnetic impurities, the effect is more dramatic. For instance, the 
transition temperature T c is reduced from its clean limit value T c0 , according to the well 
known pair-breaking formula: 

where ip(x) is the digamma function and p c = (2tttT c )^ 1 is a measure of the strength of the 
scattering and r _1 is the scattering rateii. 

By contrast in the case of superconductors whose Cooper pairs are of exotic 'p-wave' 
or 'c?-wave' character even simple potential scattering, which does not break time reversal 
symmetry, causes pair breaking!. This fact was noted already in the early contributions to 
the field0, but has become a subject of intense scrutiny only recentlyEH3. Of particular 
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interest are two dimensional models featuring 'd-wave' pairing as these may be relevant 
to experiments on high T c superconductors. Notably, for cuprates many experiments have 
explored the variation of T c , the density of states and other properties as a function of 
Ni and Zn substitutions on the copper sites@&@ or irradiation damage&"0. Although a 
wide variety of theoretical ideasS~il and methodsiHil have been applied to interpret the 
experiments a comprehensive picture of the role of disorder is far from complete. On a more 
formal level, an intriguing problem arises from the observation of Gorkov and Kalug in§ 
that the scattering in models where the order parameter has a line of zeros on the Fermi 
surface is highly singular and this may be a manifestation of interesting new physics. Indeed 
in 2-d, Nersesyan et alM predict that the quasi particle density of states N(E) approaches 
zero, even in the disordered state, as power low, ~ \E\ a , with positive exponent a, instead 
of going to a finite value 

N(0) ~ -e~7 (2) 
7 

where 7 measures strength of the interaction, as was found by Gorkov and Kalug in§. An- 
other interesting and controversial issue is the relative importance of the self-consistent Born 
approximation (SCBA) and resonant scattering in the unitarity limitffl. Our aim here is 
to explore the subject systematically on the bases of explicit calculations, albeit for a simple, 
extended Hubbard model with attractive interactions and site diagonal randomness only. 

In short we will examine the problem of disordered unconventional superconductors 
making use of the coherent potential approximation (CPA). The CPA is the most reliable 
approximation developed for the theory of electronic structure of random metallic alloys in 
the normal statefH'0. Notably it has been shown to be exact in both the weak and the 
strong scattering limits, and applicable to systems with low as well as high concentration of 
impurities. Significantly, the CPA reduces to the self-consistent Born approximation (SCBA) 
for weak scattering impurities, and agrees with the self-consistent T-matrix approximation 
(SCTA) results for strongly scattering impurities of low concentrations. Indeed it remains a 
good approximation in the unitarity limit of resonant scattering011. Finally, on account of 



the fact that it becomes exact as the number of nearest neighbors goes to infinity, the CPA 
is often referred to as a mean field theory of disorder!!. 

Given these desirable features it is clearly worthwhile to explore the consequences of the 
CPA for disordered superconductors. For the case of conventional s-wave pairing this has 
already been done, generating many useful results00. Apart from our earlier brief report@ 



and the limited discussion in Ref. [45], the case of superconductors with Cooper pairs of d 
symmetry will be treated here, within CPA, for the first time. 

We will describe fully our numerical method, demonstrate that in various limits our 
formalism reproduces many of the well known results for disordered superconductors, and 
examine in detail the phase diagram of the local and non-local attractive two-dimensional 
Hubbard models. In particular we study the variations of T c with impurity scattering 
strength and with impurity concentration for the case of local s-wave pairing as well as 
non-local (extended) s-wave and d-wave pairing. We also contrast the cases for a binary 
alloy, A-B type, disorder with the case of uniformly distributed scattering potentials on each 
site. Finally, we investigate the DOS, N(E), at low energies and its consequences for mea- 
surements of the specific heat, comparing our results with those of Kalugin and GorkovS, 
Eq. and of Nereseyan, Tsvelik and Wenge 

II. INCORPORATING THE CPA INTO THE BOGOLIUBOV DE GENNES 

EQUATION 

Our starting point is the single band Hubbard model with an attractive extended inter- 
action which is described by the following Hamiltonian 

H = 'u r in<> + \Y, UijTknj ~ 1^0 ~ ei)fk (3) 

ija ij i 

where c\ a and Cj CT are, respectively, the usual creation and annihilation operators for electrons 
on site i with spin <r, and the local charge operator is hi = h^ + h^ with fi ia = C' —iac^. The 
chemical potential is /i, Uj are the hopping integrals (for i ^ j) and e« is the local site energy. 
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The interaction term, U — ij, can be either be a local attractive interaction (Uu < 0) giving 
rise to s-wave pairing, or a non-local attractive interaction (U^ < for i ^ j) giving rise 
to <i-wave or extended s-wave pairing. Disorder is introduced into the problem by allowing 
the site energies, e$, to vary randomly from site to site. 

Starting from Eq. ([3]) we apply the Hartree-Fock- GorkovSS approximation, which 
suits in the following Bogoliubov de Gennes equation: 



re- 



E 

i 



(lUJ n —Si + jJL)8u + U 



A; 



{lUJ n + Si- fJ,)5a - U 



I 



\ (i o\ 
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IJ 



■ (4) 



G n (l,j;iu n ) Gx 2 (l,j;iu v 
j \G 2 i(l,j;iuj n ) G 2 2(l,j;iu n ) j 
for the Greens function matrix G(i,j; iuo n ) at the Mastusbara frequency, huo n = (2n+l)irkBT. 
For computational convenience we shall take the hopping integral tij to be non zero only 
when the sites i and j are nearest neighbors. The mean field pairing potential Ay can 
either be local (i = j) or nearest neighbor non-local, see Fig. 1, and couples particle and 
hole amplitudes between sites i and j. Of course, the above equations are completed by the 
self-consistency condition that 



^i j = \U ij \-Y,^ nr} G 12 {i,j;iuj n ) 

P n 



(5) 



where rj is a positive infinitesimal. To simplify matters we have assumed that the normal 
Hartree and exchange terms can be absorbed into the definitions of the chemical potential, 
/i, or the hopping integrals ty. As usual, Eqs. (fD and (|5|) are to be solved subject to the 
requirement on the chemical potential that 

2 



P n 



(6) 



where n is the average number of electrons per unit cell. Clearly, the Greens function matrix 
G(i,j; iu n ) determined by the above equations depends on the set of site energies {si}. Our 
task is to find the configurationally averaged Greens function matrix (G(i,j;tu n )). Evi- 
dently, this is made much easier if we assume that the pairing potential does not fluctuate 
from configuration to configuration. As was argued by Gyorffy et a/.0 this is a good ap- 
proximation when the T=0 coherence length £ is large. Thus our results will have to be 
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treated with appropriate care when applied to superconductors with short coherence length 
such as superconducting cuprates. 

Let us now proceed to deploy the CPA strategy for calculating the averaged Greens 
function matrix (G(i,j;iu n )) subject to the self consistency conditions: 



1 



^ij = \U ij \-J2e^(G 12 (iJ;ico n )), 

P n 



(7) 



n 



-$>^(G n (M;^n)>. 

P n 



(8) 



The first move in deriving the fundamental equations of the coherent potential approxi- 
mation is to define a coherent medium Greens function matrix G c (i,j; iu) n ) by 



E 

i 



G c (l,j\iuj n ) = 5i 



(9) 



^A u {%UJ n - fM - £ 2 2(«0>n))£»J J 

As will be clear latter G c (i,j; iu n ) = (G(i,j; iu) n )) and hence £n(m; n ) and £ 22 (2u; n ) are 
the diagonal components of the usual self-energy. Note that we did not introduce any off 
diagonal self-energies such as £i 2 (zu; n ) and S 2 i(«tu n ) because for the single site perturbations 
of our model they are zero. The next step is to consider the scattering of the quasi-particles, 
propagating according to G c (i,j;iu n ) by the defects described by the potentials: 



1 e l 1 ( E U ( H ) 0^ 







£ 2 2 («*>») 



(10) 



where I labels one of the m different site energies we wish to consider. 

In a straightforward application of the CPA principles,^ Yt{iu n ) and therefore 
G c (i,j;zuj n ) is determined by the condition that these defects do not scatter on the av- 
erage i.e. 



J2dT l (tuj n ) =0 



(11) 



where 



T l {iu n ) = V\iuo n ) l-G c (i,j;zu; n )V 



-i 



(12) 



and 

m 

5> = 1. (13) 
1=1 

From Eqs. ( |Tiyi2|) it is now possible, in conjunction with equations (0-0), to calculate H(iu n ) 
and G c (i,j; iv n ). The numerical methodology for calculating G c (i,j; iu) n ) and X(zo; n ) closely 



follows that in Ref. [48] and is outlined in Appendix A. 

A number of recent studies of superconductors with unconventional pairing suggest 
that the consequences of disorder depend sensitively on the models used to describe the 
randomness@i3 Thus we are going to investigate two different models. The first corre- 
sponds to binary alloy disorder, where m = 2. Namely we consider two types of sites with 
site energies E\ and e 2 and concentrations of c and 1 — c respectively. The second model is 
described by a uniform distribution of site energies. Here we shall have in mind the limit 
where m — > 00 with E\ G [—4, §]. Consequently in Eq. ( |TT| ) the sum J2i becomes the integral 
§ / tfei- 

In the bimodal case, where m = 2, we can simplify Eq. ([11]) to find 

E u (zu; n ) = (2c- I) 5 - -(1- E n (iu; n ))G c n (zu n )(- S ~ - S n (^ n )) (14) 
where \ei — e 2 \ =6. 

In the second, uniform distribution, case we transform the sum in equation (]TT|) into an 
integral so that 

1 ~ 

-J^T\iu n )d£i = 0. (15) 
After some straightforward algebra this leads to 

^'■^-d^+Lh^)- (16) 

Thus our CPA calculations will consist of solving numerically either Eq. (|I4]) for the 
bimodal distribution of the site energies, or Eq. ([IB]) for the case of uniform distribution, to 
determine the self-energies En(«j n ) and E 2 2(z^n)- 
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III. PAIR BREAKING FORMULA IN CPA 



We now relate the CPA formulae derived above to the usual results of disordered su- 
perconductors, corresponding to the well known pair breaking formula Eq. ([!]). As is 
well knownB the pair breaking formula was first derived for magnetic impurities in s-wave 
superconductors! but it also applies in many other interesting circumstances such as our 
present concern, namely the case of non-magnetic impurities in <i-wave superconductors0. 

To derive it within the CPA let us start with the gap equation 

A k = JjJ2 U % X £ GIM ^n)e^. (17) 

As a motivation for our argument we recall the method of Abrikosov and Gorkovi for 
solving the gap equation at T c for a clean superconductor. In that case, to find T c we linearize 



the analogue of Eq. ([17|) by approximating the off diagonal Greens function G\ 2 as follows: 



G\M zu n ) - ^ (18) 

where eg+ /i, and for our tight binding model with a square lattice = — 2t(cos(g 1 .) + 
cos(q y )). Then, we note that the kernel of the linear integral equation for A(k) is a four 
term degenerate kernel: 

U(k - q) = |E/|(f7£%+7jfc7?+ 2 sin ^ sin ^ + 2 sin k y sin q y ), (19) 

where r]^ = 2(cos(/c x ) — cos(k y )) and 7^ = 2(cos(k x ) + cos(k y )). Consequently, the general 
A(A;) will be a linear superposition of jt, sinfc^ and sink y . However, when the internal 
symmetry of the singlet Cooper pair is pure <i-wave we may take A(/c) to be of the form 

A g = A vm . (20) 

Then the condition for non-zero order parameter becomes 

l = ^E!^E^. (2!) 

Let us now define a d-wave weighted density of states 



1 



W = ^Ey^-^) (22) 



N ^ 4 

and write the above condition, which determines the transition temperature T c0 , as 

/oo 2 
deN^E)^ £ -j— -, (23) 

where u; n = nT c0 (2n +1). In the above equation the integral and the sum are divergent, so 
we need to introduce a cut-off, uj c n . In the usual way we assume the density of states N d (E) is 
slowly varying up to the cut off energy, so we will make the approximation Nd(E) = Nd(0). 
Then, considering that 



we can write 



1 = |[/|7V d (0)27rT c0 —■ (25) 



w„>0 Un 



and hence rewrite Eq. (|23|) as 

wo=*G + ^)-*G)" h (^)- (26) 

This is the BCS result for the superconducting transition temperature in the case of c?-wave 
pairingS. It differs from the conventional result only in that the d-projected density of 
states N d (0) has replaced the usual full density of states N(0). 

Let us now return to disordered superconductors and examine how the above well known 
argument is modified when the randomness is dealt with within the CPA. Using Eq. (||) it 
can be easily seen that instead of Eq. fll8|) we should use 



G c 12 (q;iuj n ) = — ^ — — --. (27) 

to linearize Eq. ([17]) at T c . Thus, noting that 

S 2 2(^n) = -En(-«o; n ) (28) 

the condition which determines T c can be written as 
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/oo 9 
dEN d (E)T c £ , - ( rr. (29) 

Now, at this point we need to know the form of Sn(zo; n ) to progress any further. As a 
first approximation we assume that the most important component to the self-energy is the 
component at the Fermi energy E = Ep = \i. Later on we will test the accuracy of this 
approximation by examining our numerical results for Yin(iuj n ). For now, however, let us 
proceed by taking 



E n (iuj n ) = i\E \sgn(u n ). (30) 



Evidently this leads to 



l = |tf|/ dEN d (E)T c Y — — — -— -. (31) 

J-oo Jr£ [iuj n - E + i\E \sgn(uj n ))(iuj n + E + i\E \sgn(uj n )) 

Again taking N d (E) outside of the integration as iV d (0) and performing the integration over 
E we find 

1 = |[/|iV rf (0)27rT c £ 1 (32) 

oj n >0 U n + 1^0 1 

where again the sum is cut off, as in the clean limit, by u^. If we now add and subtract the 
terms corresponding to S = (the clean case) we find 

1 w " 1 w " / 1 1 \ 

= 2ttT c Y — + 2vrT c Y —— - — . (33) 

Clearly the term on the LHS of Eq. ( p3|) can be replaced by In (7 2 ^t ) 011 accoun t 

of Eq. (p6|). With the same accuracy the first sum on the RHS of equation (|33| ) equals 
In (j-^A and the second sum is convergent. Hence the cutoff uj c n can be extended to 
infinity. As has been noted frequently before this infinite sum can be readily performedll 
and we find 

^ = 1> (I) ~ 4 ( I + Pc) ■ (34) 



T c q J V 2 / V 2 

where 
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Eqs. (|34"D and (|35|) are the central results of this section. Reassuringly, whilst equation 
([34]) is the standard pair-breaking formula,^ Eq. (|35| ) is a very natural, but novel, CPA 
expression for the pair breaking parameter p c . Recall that our derivation of the above 
result from CPA involved the approximation: En(zu; n ) = E . To test the validity of this 
approximation we wish to compare exact CPA numerical results with the predictions of the 
analytical expression: Eqs. (0) and fl35|) . Using numerical solutions of the CPA equation, to 
be discussed latter, Fig. 2 plots the pair breaking strength p c vs. S/t, the disorder strength 
for the binary alloy type disorder . To find pair breaking parameter p c we calculate T c for 
each disorder strength, S/t, and inverted Eq. ( [34]) to obtain an effective p c . The exact CPA 
p c can then be compared to the solid line in Fig. 2 where we have taken our numerically 
calculated values for E and directly calculated p c , via Eq. fl3~5|). Finally the dashed line in 
Fig. 2 corresponds to p c obtained using the Self- Consistent Born Approximation (SCBA). 
Evidently, the self energy at the Fermi energy, E — p = 0, E , gives a good description of the 
pair breaking parameter p c via Eq. (|34[). Also it is clear that, as expected, the Self-consistent 
Born Approximation Eq = - = tiS 2 N(0) only works well in the weak scattering limit. 



IV. ANALYTICAL FEATURES AND PREDICTIONS OF CPA EQUATIONS 

In this section we examine various analytically accessible limits of the CPA formalism 
described above. Firstly, we demonstrate that Anderson's theorem is obeyed for s-wave 
superconductors and the CPA equations are consistent with the results of Abrikosov and 
Gorkovi. Secondly we show that for ci-wave superconductors the quasi particle density 
of states at the Fermi energy, N(0), is non-zero in the presence of non magnetic disorder 
scattering and is consistent with the Gorkov-Kalugin formula, Eq. (Q). 
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A. The Anderson's Theorem in the Coherent Potential Approximation 



Formally, the CPA Eqs. (|TTD ,( |T2|) can be written in terms of renormalized Matsubara 
frequencies, u n , pairing parameter, and particle energies £r. This quantities are defined 
as follows 

u n = u n [l (36) 

V u n ) 

h = A £ (37) 

i% = ~ V ~ Re£n(2u; n ), (38) 

consequently 

Gn («*;„) = ^ E ^"tft - (39) 

k y ' ^k k 

and for the alloy type disorder with c = 0.5 and e = ±5/2, the self energy En(zo; ri ) which 
renormalizes u n , A^ and £r, and is defined by Eq. (|14D can be written as 

1 + Gfj^jEn^n) 

The alternative expression for En(?o; n ) in the case of a uniform distribution of local poten- 
tials, -5/2 < e l < 5/2 is given in Eq. (|I§. 

Note that in the case of a non isotropic, c?-wave gap, Ag, is not renormalized if the 
disorder is diagonal both in site and Nambu space. This is different from the case of local 
s-wave pairing where, as can be readily shown, A^ is renormalized by the same factor as £u n 
in Eq. (|36D . Thus for conventional superconductors, in contradiction to Eqs. (|36|) - fl39"D , we 
find that the CPA yields 

- - 1 < 41 > 

u n A 

in agreement with Born approximation or Abrikosov-Gorkov theory!. As is widely 
appreciated&0 the above equation implies the Anderson's Theorem in s-wave supercon- 
ductors. By contrast in the c?-wave case represented by Eqs. (§§), Eq. ([!]]) does not hold 
and hence there is no Anderson theorem. 
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Finally in concluding this section we would like to stress that Eqs. (|36|) - (|(]) represent 
strictly pure ci-wave result. Even if we stick to the singlet case a more general solution of the 
CPA equation will imply a renormalization of Ag to A^. A good example of such situation 
is a case where the symmetry of A^ is of extended s-wave symmetry s* oc (cos(k x ) + cos(k y )) 
type. We shall encounter this interesting circumstance later on in this paper. 



B. Density of States iV(0) in d wave superconductors 

Moving on and returning to the ci-wave case, we observe that the form of Eqs. fl36|) - 
([40|) are the same as were found by Larkin.0 Thus again the CPA reproduces the expected 
general form of the gap and frequency renormalizations, but with an improved description 
of the disorder. The most prominent feature of a conventional superconductor is vanishing 
of the quasi particle density of states N(e) for energies e measured from the Fermi energy 
Efi less then A. In the case of clean, ci-wave superconductors the line of zeros of A(k) on 
the Fermi surface leads to finite N(e) for all e except e = 0. In fact as is well known@ N(e) 
approaches zero linearly in e. In the present section we shall investigate what happens to 
N(e) in the presence of disorder. 

As it turns out for a given gap parameter A& = At^ and in the limit of small disorder 
5^0 the CPA equations can be solved analytically. To affect the solution note that in 



Eq. fl39|) the major part of the summation is coming from the four singular points in the 



Brillouin zone where the denominator vanishes: 

iu n 2 - H - A 2 , = (42) 

Linearizing around these points and performing the summation over k analytically, we find 
that 



(4A) 2 + (ImS 11 (0)) i 



(ImS n (0)) 2 

where a = 2tAn. Clearly in the limit |ImSn(0)| << 4A this leads to 



(43) 



14 



gu(0) ^ Im£ll(0) ln 



4A 



(44) 



ImSn(O) 

Moreover, when Hn(jw n ) is small compared to the band width we can rewrite Eq. (|40|). 
using the Self-Consistent Born Approximation as 

5 2 

E u (iu n ) = —G c n (iu n ). (45) 



On substituting this result into Eq. (|44|) the later becomes an equation for (5^(0) which 
determines the density of quasi particle state N(0) via the formula A(0) = ^G^iiUn = 0). 
Indeed we find 

» T /„N 4At 87rAt 

and 

- ImEn (0) « 4 Ate" ^ . (47) 

These formulae agree with the results of Kalugin and Gorko^ and Haas et al£5 and 
have been verified numerically. For example Figs. 3 and 4 compare N(0) as calculated 
from Eq. ([46]) and calculated using completely self-consistent CPA. As one can see there is 
good agreement between the two results and thus we conclude that at e = the density of 
states, in CPA, becomes finite when an arbitrary small amount of disorder is introduced. 
However, one should note that at very low levels of disorder our numerical results tend to 
a constant , whereas the analytical results tend to zero, this comes from the small complex 
component we have added to the energy to enable us to evaluate the Greens function. For 
a more detailed description of what is happening to iV(0) and — ImEn(O), in the presence 
of disorder, see Sec. VII of this paper, where we have analyzed their properties more closely 
and report, extensively, further numerical results. 



V. LOCAL QUASI-PARTICLE DENSITY OF STATES CALCULATIONS 



In this section we present results for the single site local quasi particle density of states 
calculations. As mentioned before we consider two types of disorder, (i) binary alloys where 
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we have two on-site potentials randomly distributed throughout the lattice and (ii) a uni- 
form distribution of random on-site potentials. For both types of disorder we have solved, 
numerically, Eq. ([5p in conjunction with its appropriate self-consistency conditions, for the 
order parameter, A, the average number density, n, and the self-energy, £n(Kj n ). 

The first situation we consider is a binary alloy of random on site energies E\ and e<i 
with equal concentration c = 1 — c = 0.5. The parameter we have chosen to use to describe 
the strength of the disorder is 5 = E\ — Figure 5(a) shows the density of states in the 
normal state. The Van Hove singularity characteristic of a tight-binding model with nearest 
neighbor hopping on a square lattice is clearly visible for small disorder (5 = 0.6t) in the 
middle of the band. As Fig. 5(a) shows, for more strongly disordered alloys this Van Hove 
peak is split into two peaks with some additional smearing. In the limit of 5 being very large 
we get, as one would expect, band splitting of states associated with e\ and e 2 respectively. 
On the other hand the disorder with uniform distribution of site energies, Fig. 6(a), gives 
only the smearing and flattening of density of states with, eventually, a complete flattening 
of the Van Hove peak. 

The imaginary part of the self-energies for these two types of disorder are shown in 
Figs. 5(b) and 6(b) respectively. Evidently at low levels of disorder, they reflect the density 
of states of the clean system, and hence are consistent with the Born Approximation via. 
—lmEii(E) oc N(E). For larger disorder strengths ImS n (£) is more flat for the uniform 
distribution, whereas for the binary alloy the self-energy is peaked not at the Van Hove 
singularity but at S\ and e 2 - 

Before turning to the problem of disordered <i-wave superconductors for reference we 
have studied, briefly, the s-wave case. In short, we have introduced a site diagonal, local, 
attractive interaction with strength U, into the above model. As expected such interaction 
leads to conventional s-wave pairing. When implementing the CPA we have assumed that 
in this case the pairing order parameter A, which is now site diagonal does not vary from 
site to site. This is consistent with the Anderson theorem which was shown to be adequate 
for systems with a long coherence length.lll Fig. 7(a) shows the s-wave quasi-particle density 
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of states as calculated for various values of the binary alloy disorder strength 5. The results 
confirm that, for s-wave pairing, the gap is absolute and, whilst the edges may move, 
they remain well defined as required by Anderson's Theorem^. In Fig. 7(b) we have also 
plotted the self-consistent self-energy as a function of the quasi-particle energy, for the same 
disorder strengths used to obtain Fig 7(a). From this graph we can see that -ImEn(E) is 
zero inside the gap and hence the normal disorder does not act as a pair breaker, for s-wave 
superconductivity. 

The quasi-particle density of states for <i-wave superconductors is dramatically different 
from the above BCS spectrum in the s-wave case even in the clean limit. As it is well 
known it has the characteristic v-like dip shown in Fig. 8(a), for e near the chemical 
potential (/i = 0). Upon introducing binary alloy disorder into the <i-wave system one 
would naturally expect similarly different behavior from that described above and in Fig. 
7. As is clear from the results reported in Fig. 8 this is indeed the case. Strikingly, as the 
analytic results of the previous section suggested, N(0) becomes non-zero for the slightest 
disorder. This implies gapless superconductivity in contrast to the gaped one in the s-wave 
case. 

In Fig. 8(b) we have plotted the self-energies, corresponding to the quasi-particle density 
of states results presented in Fig. 8(a). This shows a complex evolution with disorder, for 
small 5 the imaginary part of the self-energy, reflects the pure ci-wave density of states 

(as expected in the SCBA limit) and hence Im[E(0)] ~ 0. However increasing S leads to a 
finite E(0), with a cusp like minimum in Im[£(i?)] at E = 0. Increasing the disorder even 
further, to 5 = 2.8t, the <i-wave pairing is completely destroyed, and ImE(E) reverts to the 
normal system self-energy. In this case — ImE(E) is a maximum at E = 0, since the Fermi 
energy, E F = fi, was set exactly to 0. 



17 



VI. CRITICAL TEMPERATURE CALCULATIONS 



To analyze the effect of the disorder upon T c we have solved the gap equation together 
with the CPA equations. Again for the sake of comparison we have developed the analogous 
theory for the conventional s-wave superconductors based on the site diagonal particle- 
particle interaction, of strength Uq, mentioned briefly in a previous sections. In this case 
neglecting the spatial fluctuations of A i; and linearizing the gap equation for the configura- 
tionally averaged single site pairing potential A we find the condition for the temperature 
T c , below which A 7^ 0, to be 

where Ni(E) is the averaged density of states in the normal state at energy E, as calculated 
within the CPA, and (3 C = -^F- 



The solution of Eq. ( f48|) for T c as a function of band filling and strengths of alloyed 
disorder, 5 is shown in Fig. 9. From this figure we can see that the superconducting state 
exists for all band filling and the strength of disorder does little to suppress T c . Again this 
is consistent with Anderson's theorem. 

Now let us turn to the case where the interaction is non-local and the pairing potential Ajj 
connects nearest neighbors sites. Near T c , where the gap equation is linear, the solutions can 
be labeled by their symmetries. Indeed we find two separate conditions on the temperature 
T c for the instability of the normal state to d and extended s symmetry breaking. In the 
first case T c d is determined by 

1- lUl 1 V r JE Im ( E) \ tanh ( % E \ (40) 

1 --^n\I-o, dE Im \2E-^{E)-^E) ) tanh ) (49) 

and in the second extended s-wave case, T c s is given by 



u% r dE [m ( fffQ r \ tanh (S£ 

4tt N Y J-00 \2E - Z n (E) - S 22 (E) J \ 2 / 



(50) 



2E-Z 11 (E)-X 22 (E) i 
where rjj: and 7^ refer to s-wave and <i-wave like 'harmonics' defined previously following 
Eq. m 



In Fig. 10 the critical temperature for both <i-wave and extended s-wave pairing is 
shown as a function of band filling, n, for various strengths of alloyed disorder, 5. Full lines 
correspond to d-wave T c while dashed ones to extended s-wave solutions. In the clean limit 
(5 = 0) we can see that the extended s-wave solution exists mainly at the band edges and 
the (i-wave solution is largely confined to the central portion of the band. For the interaction 
strength we have chosen, the two solutions cross at n ~ 0.38. Evidently for n < 0.38 the 
superconducting instability is at T c s whilst n > 0.38 the transition temperature is T c d As we 
increase the strength of the alloyed disorder (S/t = 0,1.2,2.0,2.5,2.6,2.7,2.8,3.0 from top 
to bottom curve) we see that both the T c for the extended s-wave and <i-wave solutions is 
reduced and for particularly strong disorder (5 = 2.7t and S = S.Ot) the maximum in the d- 
wave T c is no longer at n — 1. This is connected with the splitting of Van Hove singularities 
visible in Fig. 5(a), i.e. the maximum in the normal-state density of states corresponds to 
the maximum in the d-wave T c d . 

On the other hand for a uniform distribution of random site energies where there is 
no splitting of Van Hove singularities, and hence, we see a different tendency (Fig. 11, 
S/t = 0, 1.0,2.0,3.0,4.0,5.0,6.0 from top to bottom curve). The maximum value of T c d is 
located at n = 1 no matter what the strength of disorder. For extended s-wave solutions 
and for relatively small disorder in the of d-wave case the decreasing tendency in T£' d with 
growing disorder are very similar for the two types of disorder. However above a certain 
strength in either case we note that there is no crossing between the extended s-wave and 
o?-wave solutions. This fact may be interpreted as the sign that the s and d sates can not 
coexist or mix. 

To compare the effect the two types of disorder have upon T c we have, in Fig. 12, 
calculated T c d vs. S/t at half filling for both binary alloy (solid line) and uniform (dashed 
line) disorder. From this graph we can see that the two types of disorder effect T c in almost 
exactly the same manner. If we rescale the graphs such that the standard deviation in the 
disorder strength is the same, for uniform disorder a = and for binary disorder a — |, 
then the two curves lie very close to one another. This prompts the conclusion, which may 
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however be premature ,0 that T c is not very sensitive to the details of the fluctuation in the 
site energy. 

Other interesting point to investigate concerns the relative robustness of T c against degra- 
dation by disorder in the cases of conventional s and <i-wave pairing.0 Fig. 13 shows the 
results for T c for inter-site c?~wave and on site s-wave superconductors versus alloyed disor- 
der strength 5, where n = 0.6. Here one can easily recognize a typical difference between 
these two superconducting states. Clearly in case of a <i-wave superconducting state disor- 
der acts as an effective pair breaker while for s-wave superconductors it decreases T c only 
slowly or not at all. An interesting physical consequence of this effect is that if both Ua 
and Uij (i ^ j) are attractive increasing disorder could lead to a transition from d-wave to 
s-wave pairing, as suggested by AbrikosovH. 

Evidently, the binary alloy disorder has certain features not shared by the uniform dis- 
tribution of site energies model. One of these is the possibility of band splitting. For a very 
large 5 parameter the band is split into two sub-bands, as shown in Fig. 14. The on site 
s-wave superconducting state is created within one of sub-bands. For disorder parameters: 
c = 0.5 and t\ = —t^ = 5/2 we have not found any similar solution for c/-wave, but for 
c=0.25 and n=0.5 a single sub-band rf-wave solution is found. In Fig. 15 we show the very 
different behavior of T c for a small and a large non-local interaction parameter, U^. For 
large enough U^ the ci-wave T c does not go to zero when the band splits. The quasi-particle 
density of states corresponding to such a possibility is presented in Fig. 16 and, as in Fig. 14, 
we see the superconducting gap, which is now d-wave, has formed in one of the sub-bands. 

To complete the discussion Figs. 17, 18 and 19 show the critical temperature plotted 
versus concentration c for various strength of binary alloy disorder 5 with band fillings of 
n = 1 (Fig. 17), n = 0.7 (Fig. 18) and n = 0.1 (Fig. 19). Figs. 17 and 18 correspond 
to (i-wave while Fig. 19 correspond to inter-site extended s-wave pairing. Again for large 
enough alloy disorder pair breaking phenomena take place. If disorder is introduced by a 
fixed potential 8, but varying the concentration of scatters, c << 1, from we see that the 
results are similar to the case of a fixed c and increasing 5 from 0. This shows that CPA 
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can work equally well, and indeed asymptotically exact, in both the Born (small 5) and 
T-matrix (small c) and resonant scattering regimes. 

VII. N(E) AND Y, U (E) AS E -> AND THE LOW TEMPERATURE SPECIFIC 

HEAT 

As mentioned earlier the behavior of the quasiparticle density of states N(E) and the 
imaginary part of the self energy ImSn(E) near E = is of general conceptual significance. 
For instance, the power law behavior of N(E) oc \E\ a , for ci-wave superconductors give 
rise to power law dependence with temperature of many thermodynamic quantities, such 
as the specific heat c v (T), instead of the exponential cut-off characteristic of a gap in the 
quasiparticle spectrum.! Naturally, dramatic changes in the low energy behavior of N(e) 
and ImEn(e) as disorder is added to the problem is also of general interest and, as it turns 
out, of lively controversy.llHl3E3'il In this section we wish to examine those predictions of 
the CPA calculations which are relevant to these issues in particular. Using the methods 
outlined in the previous sections we calculate numerically the quasi-particle local density of 
states, N(e), and the self-energy, ImS n (e), and then investigate how these two quantities 
change with both temperature and disorder. 

For simplicity we have studied the half filled band (fi = 0, n = 1). This is the band 
filling for which T c , in the <i-wave case, is a maximum. We will examine the effects of alloyed 
disorder on the system and show that as we increase the amount of scattering the specific 
heat vs. temperature relation changes from T 2 , in the clean limit, to T in the limit of 
strong scattering. As is widely recognized! this behavior is the consequence of N(E) ~ \E\ 
changing to N(E) = constant and is consistent with experiments.! 

To get an impression of the form of N(E) and — ImSn(E) in the region of the chemical 
potential, /x, we fit them to the function a + b\E\ c in the energy range —AE < E < AE, 
where AE is small compared to the gap. Using the coefficients a, b and c to fit the curves 
N(E) and —ImHu(E) in the region of \i gives us a tool to analyze their functional form. 
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For example a tells us if the curves are finite at E = 0, b controls the gradient of the curve 
and c controls the curvature. 

In Figs. 20(a-c) we have plotted these coefficients, for N(E), as a function of T, for 
different values of 5. We also included in Fig. 20(d) plots of the magnitude of the d- 
wave superconducting order parameter |A| vs. T, for the same values of alloyed disorder 
strength, 5, to indicate the temperature where the order parameter goes to zero. Regarding 
these curves as a brief summary of what the CPA predicts about the low energy behavior 
of N(E) and we now comment on their implications. 

Firstly we note that in Fig. 20(a) the parameter a(T) tends to a finite limit as T — > and 
this limit increases more and more rapidly as the disorder described by 5, increases. Hence 
N(0) is finite in agreement with our earlier discussion in Sec. IV where we derived for low 
scattering and low temperatures the dependence of N(0) on S (see Eq. (f46|)). It lends credit 
to the general consistency of our results that the a(T) = N(0, T) curves rise to their normal 
state values as T —>■ T c . Interestingly, as can be seen in Fig. 20(b) the gradient of N(E), 
namely b, changes dramatically only near T c . In fact for the case where S is small (squares) 
the gradient even changes sign. This means that we go from having a gap in the density 
of states, below T c to having a Van Hove singularity above T c . The curvature, represented 
by c and plotted in Fig. 20(c), shows that for low disorder (squares) as we increase the 
temperature we go from a curvature of c < 1, ie. a cusp, to a curvature c > 1. Finally when 
the critical temperature is reached for each of the four different disorder types c goes to 2. 

In the same manner, we now examine the corresponding coefficients, a', b' and d for 
— ImSn(£/). In Figs. 21(a-c) we have plotted the calculated coefficients for -ImSn(£), at 
different values of 6 as a function of the temperature. In Fig. 21(a) we can see that at low 
temperatures — ImSn(0) increases with disorder, again this agrees with our results in Sec. 
IV where an analytical form for the dependence of — ImSn(0) upon 5 was derived at T = 0, 
see Eq. P7|). As we increase the temperature we can see the self-energy at the chemical 
potential also rises until the normal state value is reached at the critical temperature. In 
Fig. 21(b) we can see how the gradient changes from being large for large S and small for 
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small 8 at low T, and at temperatures greater than T c to be large and negative for large 
S and small and negative for small 5. Finally Fig. 21(c) shows how, c, the curvature goes 
from being almost linear, c = 1, at low temperatures, regardless of 5, to c = 2 above T c . 

To illustrate the consequences of these results for the low temperature thermodynamic 
properties we have calculated the effect of disorder upon the specific heat. Using the above 
temperature dependent coefficients for the density of states near e ~ we calculated the 
limiting behavior of the specific heat as T — > 0. The results are shown in Fig. 22 for different 
strengths of disorder 5. In the case where 5 is small (squares) we can see that the specific 
heat has a T 2 dependence (the plotted continuous solid line with no points) and for large S 
the dependence upon temperature is linear, as expected.! 



VIII. CONCLUSIONS 



We have compared and contrasted the effect of disorder on conventional s- and rf-wave 
superconductors on the bases of an extended, negative U Hubbard model and a mean field, 
CPA, treatment of disorder. On the one hand we have derived many of the well known 
results, such as the pair breaking formula in Eq. ( pM| ) or that for the quasi-particle density 
of states iV(0), Eq. (f46|) . On the other hand we have solved the Gorkov-CPA equations 
Eqs. fll4]), ( |17D numerically and surveyed the salient features of their consequences by 
explicit calculations. The use of CPA in describing disordered <i-wave superconductors is 
an advance in this very active because it allows us to avoid the usually delicate 

choice between methods, sets of diagrams, designed to deal with either weakly or strongly 
scattering perturbations. The CPA treats both kinds of problems equally accurately and it 
is known to provide a very credible interpolation between the twc@. As an example where 
above feature of CPA may have a crucial role to play we recall the use of the resonant 
scatterer model in interpreting experimental data both on the cuprates&S and some Heavy 
Fermion systems.0'0 In short, note that in the impurity, c — > 0, limit the self energy for the 
Greens function describing an electron moving on a lattice, is given by 
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E(E) = cT(0,0;E), (51) 

where the T-matrix is defined, in terms of the impurity 'potential' Vi mp , a real local shift in 
the site energy, as 

Evidently T(0, 0; e) is complex number with an amplitude: 



a(E) = ^(1 - V mp ReG(0, 0; E)) 2 + (V imp lmG(0, 0; E)) 2 

and a phase, the phase shift <f)(E), given by 

V, mr ImG(0,0;E) 
tan ^ )= l^V„„ p Re G (0,0; E ) - (M) 

Now, observe that while for weak scatterers 

Im£ B (£) = -C7T | V 4mp | 2 iV(£) (54) 



for a resonant scatterer in the unitary limit, defined by 4>{E) 



2 ' 



hnE*{E) = -c^, (55) 

where N(E) = — -ImG(0, 0; E). It is the above dramatic difference in the dependence of 
ImS s and ImE^ on the density of states N(E) that the cited authors rely on in interpreting 
the relevant experimental data. Evidently, since the individual scattering events described 
by the local T-matrices are always treated exactly in the CPA, the CPA describes weak 
scatterers and resonant scatterers equally well. Moreover, since it is a reliable approximation 
not only in the impurity limit, c ~ 0, but also for arbitrary concentrations it deals with 
resonance scatterers even when Eq. (|55D no longer holds. Thus CPA should be the preferred 
treatment for models with strong, even resonant scattering. 

Although we have not specifically concentrated on this aspect of the method, the principle 
feature of Eqs. (|54]) and (|55|) , namely the dependence of T,(E) on N(E), can be seen to be 
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at work in our calculations of the previous section. To demonstrate this we have calculated 
S(0) and N(0) as functions of band filling, n, in the most interesting region n ~ 1.0, and 
compared in Fig. 23 their relationships for weak and strong scatterers at c = 0.5. Clearly, 
for weak scatterers 5/t = 1, — ImS(O) is a more or less linear function of N(0) as in Eq. 
(0), while for strong scatterers 5/t = 3, ImS(O) is inversely proportional to N(0) as in the 
resonant scatterers model described by Eq. (f)5"l). Thus we conclude that the CPA employed 
in the calculation gives a reliable account of disorder in both the weak and strong scattering 
regimes. 

Having listed the above desirable properties of the CPA we hasten to emphasize that it is 
a 'mean-field' theory of disorder and, hence, does not describe such interesting phenomena 
as localization!! even in the normal state.0 Consequently, our result that N(0) is finite 
for the smallest amount of disorder can not be taken as evidence against the conclusion 
that N(E) ~ \E\ a of Nersesyan et aL0 As this originates from the divergence of the 
vertex corrections in perturbation series for T,(E) we may conjecture that it has to do with 
localization effects not described by CPA. This very interesting point is in the need of further 
clarification, and indeed we shall return to it in a later publication. 



ACKNOWLEDGMENTS 



This work has been partially supported by EPSRC under grant number GR/L22454 and 
KBN grants 2P 03B 031 11 and 2P 03B 050 15. A.M. Martin would also like to thank the 
TMR network Dynamics of Nanostructures. 



APPENDIX A: THE RECURSION METHOD FOR CPA 

The recursion method@ was first used for solving the Bogoliubov de Gennes Equation 
by Annett and GoldenfeldS and Litak et a£3. Its use to implement CPA calculations was 
advocated by Julien and Maycll. Here, as in reference^ we made a combined use of these 
powerful methods to calculate the Greens function matrix 
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•b'a'> 



(Al) 



'El - H 

where the indices i and j denote sites, while a and a' represent the particle or hole degree 
of freedom on each site. We denote particle degrees of freedom by a = 1 and hole degrees 
of freedom by a = 2. For example G\ 2 (i, j; E) represents the Greens function between the 
particle degree of freedom on site i and the hole degree of freedom on site j at energy E. 

To compute the Greens functions ( |A1|) we can closely follow the method described by 
Martin and Annett.E§ Using this method we can transform the Hamiltonian to a block 
tridiagonal form 



El-H 



El — a —61 
-b\ El - ai 





-b 2 














-bl 

































V 



(A2) 



-bl El - a r 



where a n and b n are 2x2 matrices. Given this form for (ia\El — H\ja') and expressing 
the Greens function as 



G c aa ,^3;E) = (ta\(El-H)- 1 \ja') 



the Greens functions above can be evaluated as a matrix continued fraction so that 



(A3) 



G(i,j;E) = I El - a - b{ [El - ai -b\ (El - a 2 - b\ (El - a 3 



l b, 



-1 



-1 



-1 



(A4) 



where 



G c (i,j;E) 



/ G c aa (i,i;E) G° aa ,(i,j;E) ^ 



V 



(A5) 



G c a , a (j,i;E) G c a , a ,(j, r ,E) 

Within Eqs. ( |A2] ) and (|A3j) we have a formally exact representation of the Greens 
functions. However in general both the tridiagonal representation of the Hamiltonian, and 
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the matrix continued fraction ( |A2|) will be infinite. In practice one can only calculate a 



finite number of terms in the continued fraction exactly. In the terminology of the recursion 
method it is necessary to terminate the continued fraction-iUBii^ii 

If we were to calculate up to and including a n and b n and then simply set subsequent 
coefficients to zero then the Greens function would have 2n poles along the real axis. The 
density of states would then correspond to a set of 2n delta functions. In order to obtain 
accurate results it would be necessary to calculate a large number of exact levels, which 
would be expensive in terms of both computer time and memory. 

As a more efficient alternative we choose to terminate the continued fraction using the 
extrapolation method, as used previously by Litak, Miller and Gyorffy.E3 We calculate the 
values for a n and b n exactly up to the first m coefficients using the recursion method. Then, 
noting the fact that the elements of the matrices a n and b n vary in a predictable manner ,H 
we extrapolate the elements of the matrices for a further k iterations, where k is usually very 
much greater than m. This enables us to compute the various Greens functions accurately 
with relatively little computer time and memory. 

We now have a method for calculating the Greens Function G^ a ,(i, j; E). To calculate 
our system of equations self-consistently for a given site i, we need to know G\ 2 {h3i] E), 
G\ 2 {.h32] E), G\ 2 {i,jz] E) and G^ihji', E). If the interaction is non-local, else if the inter- 
action is local then we only need to calculate Gi 2 (i, i', E), see Fig. 1. To calculate n we also 
need G\ l {i ) i; E). 

In performing our self-consistency calculation we can make use of the fact that S(-E') just 
acts as a renormalization to E i— > E (see Sec. IV). This means we do not need to recalculate 
a n and b n when we calculate the self-consistent self-energy. We first calculate the order 
parameter self-consistently for a given attractive interaction, then we calculate G c (i,i;E) 
and hence the self-energy self-consistently (remembering that in this self-consistent loop we 
do not need to recalculate a n and b n , since the self-energy does not affect them). Then 
we recalculate the superconducting order parameter self-consistently, incorporating the cal- 
culated self energy and repeat this procedure until both the order parameter and the self 
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energy have reached stable self-consistent solutions. 
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FIGURES 

FIG. 1. Schematic diagram of a tight binding lattice with particle and hole degrees of freedom. 
The gap parameter Ajj couples particles on site i to holes on site j. The difference between local 
and non-local pairing is highlighted by the dashed (local pairing) and solid (non-local pairing) lines. 

FIG. 2. The effective pair-breaker p c as calculated (i) by numerically finding ^ and inverting 
Eq. (36) (squares), (ii) numerically finding ImSo and using this in Eq. (37) (solid line) and (iii) 
Using the self-consistent Born approximation to find ImSo and then evaluating Eq. (37) (dashed 
line). 

FIG. 3. A comparison of the density of states at the chemical potential vs. different strengths 
of alloyed disorder 5/t. The solid line is the analytical form derived in Eq. (48) and the dashed 
one represent are our self-consistent numerical calculations. 

FIG. 4. A comparison of the density of states at the Fermi energy vs. different strengths of 
disorder with a uniform distribution between — | and |. The solid line is the analytical form 
derived in Eq. (48) and the dashed one represent our numerical calculations. 

FIG. 5. (a) Density of states N(E) and, (b), self-energies ImE(E) for a normal state with 
various alloyed disorder strengths 5/t, calculated for a local interaction \U\ = 3.5t and n = 1. 

FIG. 6. (a) Density of states N(E) and (b) self-energies ImE(E) for a normal state with various 
uniform disorder strengths 5/t, calculated for a local interaction \U\ = 3.5i and n = 1. 

FIG. 7. (a) Density of states N(E) and (b) self-energies ImS(E) for a superconducting state 
with various alloyed disorder strengths 5/t, calculated for a local interaction \U\ = 3.5i and n = 1. 

FIG. 8. (a) Density of states N(E) and (b) self-energies ImE(.E') for a superconducting state 
with various uniform disorder strengths 5/t, calculated for a non-local interaction \U\ = 3.5t and 
n = 1. 
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FIG. 9. Critical temperature T c vs. band filling n for a local s-wave superconductor with 
various alloyed disorder strengths 8/t, calculated for a local interaction \U\ = 3.5i). 

FIG. 10. Critical temperature for ci-wave and extended s-wave order parameters as function 
of band filling for different strengths of alloyed disorder 8/t, calculated for a non-local interaction 
\U\ = 3.5t. 

FIG. 11. Critical temperature for ci-wave and extended s-wave order parameters as a function 
of band filling for different strengths of uniform disorder 8/t, calculated for a non-local interaction 
\U\ = 3.5t. 

FIG. 12. Critical temperature of the ci-wave superconducting state vs. the strength of disorder. 
The solid line is for a uniform distribution of disorder and the dashed line is for a binary alloy, 
where c = 0.5 and in both cases n = 1 and the non-local interaction \U\ = 3.5i. 

FIG. 13. Critical temperature of ci-wave (solid line) and local s-wave (dashed line) supercon- 
ducting pairing states as a function of alloyed disorder for \U\ = 3.5i, n = 0.6 and c = 0.5. 

FIG. 14. Local density of states, at zero temperature, for an alloyed disorder of 8/t = 6, a local 
interaction \U\ = 3.5i, n = 0.6 and c = 0.5. 

FIG. 15. Critical temperature of the ci-wave superconducting state as a function of alloyed 
disorder strength shown for different interaction strengths, c = 0.25 and n = 0.5. 

FIG. 16. Local particle density of states, for an alloyed disorder of 8/t = 8, a non-local inter- 
action \U\ = 3.5£, n = 0.5 and c = 0.25. 

FIG. 17. Critical temperature vs. concentration, at different alloyed disorder strengths, for a 
non-local interaction \U\ = 3.5i and n = 1.0. 

FIG. 18. Critical temperature vs. concentration, at different alloyed disorder strengths, for a 
non-local interaction \U\ = 3.5i and n = 0.7. 
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FIG. 19. Critical temperature vs. concentration, at different alloyed disorder strengths, for a 
non-local interaction \U\ = 3.5t and n = 0.1. 

FIG. 20. The coefficients a (a), b (b) and c (c) which represent the Local Particle Density of 
States in the region of the chemical potential as a function of temperature, for different alloyed 
disorder strengths S (S = 0M (diamonds), 5 = l.Ot (+), S = 2.0 (squares) and 5 = 2M (x)). In 
(d) we have also plotted the magnitude of the d-wave gap vs. temperature, for the same alloyed 
disorder strengths as in figures (a), (b) and (c). In all four figures the interaction is non-local with 
\U\ =3.5t, n= 1 and c = 0.5. 

FIG. 21. The coefficients a! (a), b' (b) and d (c) which represent the imaginary part of the 
self-energy, in the region of the Fermi energy, as a function of temperature, for different alloyed 
disorder strengths 5 (5 = O.Qt (diamonds), 5 = l.Ot (+), 5 = 2.0t (squares) and 5 = 2.Qt (x)). The 
interaction is non-local with \U\ = 3.5i, n = 1 and c = 0.5. 

FIG. 22. The specific heat as a function of temperature for different strengths of alloyed disorder 
5 (5 = 0M (squares), 5 = l.Ot (circles), 5 = 2.0t (triangles) and 5 = 2M (+)), again the interaction 
is non-local with \U\ = 3.5t, n = 1 and c = 0.5. 

FIG. 23. Relations between ImS(0) and N(0) for weak S/t = 1.0 an strong S/t = 3.0 scatterers 
in binary alloys (c = 0.5). 
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